Universal Critical Velocity for the Onset of Turbulence of Oscillatory Superfluid Flow 
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The critical velocity v c for the onset of turbulent drag of small spheres oscillating in superfluid 4 He 
is frequency dependent (uj/2h from 100 Hz to 700 Hz) and is described by v c — 2.6 i/ku, where k is 
the circulation quantum. A qua litative analysis based on a recent theory of the onset of superfluid 
turbulence gives v c « \J%KUj/j3, where /3 ~ 1 depends on the coefficients of mutual friction. This 
agrees well with the data and implies that v c is a universal critical velocity that is independent 
of geometry, size, and surface properties of the oscillating body. This is confirmed by comparing 
our data on spheres with v c obtained with other oscillating structures by other groups. Numerical 
simulations indicate somewhat larger critical velocity, above which a rapid increase in vortex length 
is observed. 

PACS numbers: 67.25.dk, 67.25.bf, 47.27.Cn 



INTRODUCTION 

A quantitative description of the critical velocity for 
the onset of turbulence caused by a macroscopic body 
moving in superfluid helium is still an unsolved problem. 
At low velocities the superfluid flowing around the body 
exhibits pure potential flow and the only drag force is due 
to scattering by thermally excited quasiparticles (rotons 
and phonons in 4 He) whose density rapidly vanishes when 
the temperature is decreased to the millikelvin regime. 
If, however, the velocity is increased a sudden onset of 
a large and nonlinear drag force is observed at a criti- 
cal velocity v c . This is due to a breakdown of potential 
flow and the onset of turbulence due to vortex produc- 
tion by the moving body. These vortices all have the 
same quantized circulation k = h/m « 10~ 7 m 2 /s (h is 
Planck's constant and m is the mass of a 4 He atom) and 
a core of atomic size. These vortices form a tangle that 
is described as "superfluid turbulence" or "quantum tur- 
bulence" . There is a large body of experimental evidence 
that v c depends on the size of the moving body [1]. While 
this is true for steady flow, the question arises whether 
the physics will be different for oscillatory flow. After 
all, for practical reasons, often an oscillating body like a 
sphere, a vibrating wire, a small tuning fork, or even a 
vibrating grid are employed in the experiments 0. 

In our present work, we first show experimentally that 
in the case of oscillatory flow generated by an oscillating 
body there is a universal critical velocity, that is given 
only by the oscillation frequency u>, the circulation quan- 
tum k, and a dimcnsionless number of order one, i.e., 
v c ~ y/ k ■ uj. This result can be made plausible simply 
on dimensional grounds Finally, an analysis based 

on a rigorous theory for the onset of superfluid turbu- 
lence in steady counterflow, when extended qualitatively 
to the case of oscillating flow in our experiments, gives 
good agreement with the data, including their very weak 



temperature dependence from millikelvin temperatures 
up to 1.9 K. 

EXPERIMENTS 

The experiments with spheres were performed at Re- 
gensburg University several years ago. Most of the re- 
sults have been published in earlier work [5j, ly, |2| , but a 
systematic analysis of the frequency dependence of the 
critical velocities was only performed recently. The data 
were taken with 2 different spheres, having a radius of 
100 /im and of 124 /xm. The smaller one was used in a 
3 He cryostat, the other one in a dilution cryostat. The 
magnetic spheres were placed inside a horizontal par- 
allel plate capacitor consisting of niobium electrodes to 
which a voltage was applied while the temperature was 
dropped below the superconducting transition tempera- 
ture of niobium, ca.10 K. At helium temperatures the 
charged spheres were levitating between the capacitor 
plates. Frozen flux in the electrodes was necessary for 
a stable vertical oscillatory motion when a small ac volt- 
age was applied at resonance. The frequency of the os- 
cillators could not be controlled quantitatively because it 
depended on the frozen flux of the particular levitation 
status. It was observed, however, that a faster cooling 
rate from 10 K to 4.2 K in general lead to higher reso- 
nance frequencies. The critical velocities were measured 
by reducing the driving ac voltage in the turbulent flow 
regime (identified by the large nonlinear turbulent drag) 
until turbulence vanished. This procedure was necessary 
because if, instead, the drive was increased from zero, it 
was found that often the critical velocity could be sub- 
stantially exceeded before potential flow eventually would 
break down 0, see Fig. Q] The origin of this hysteresis 
will be discussed in the following Chapter. 

At temperatures below 0.5 K the hysteretic behavior 
is replaced by an instability of turbulence leading to an 
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FIG. 1: (color online) Velocity amplitude as a function of 
the driving force for a 100 /im sphere oscillating at 236 Hz 
in superfluid 4 He at 1.90 K. Red circles: data taken with 
increasing drive; blue x : decreasing drive. Because of the 
strong hysteresis v c — 46 mm/s can only be determined with 
the down sweep. 



intermittent switching between potential flow and tur- 
bulent phases whose lifetimes gradually vanish when v c 
is approached by reducing the drive [5j. This instabil- 
ity has been attributed to statistical fluctuations of the 
length of vortex lines per unit volume L Q . The resulting 
v c of the various runs for both spheres are displayed in 
the Fig. [2] together with critical velocities obtained with 
other oscillating structures by other groups. 

We first discuss our own results for the spheres. The 
data in Fig. [2] were all taken below 1 K where the criti- 
cal velocities were found to be independent of tempera- 
ture, see below. The errors of the individual values for 
the spheres vary from ca. 1% to about 5%, depending 
on how closely v c was approached in the experiments. 
The exceptionally large scatter of the three data points 
near 200 Hz is by far larger than any error bar, the rea- 
son of which is not clear: only a small systematic error 
could have occurred in the calibration of the velocity scale 
(caused by a small error of few percent in the determi- 
nation of the electric charge of the oscillating sphere). 
The mostly higher resonance frequencies of the smaller 
sphere are due to its smaller mass and probably also due 
the usually faster cooling rate from 10 K used with the 
3 He cryostat. 

The other data we could find in the literature are the 
following. 1. A result from an experiment with a vibrat- 
ing wire of the Lancaster group [9( . In that experiment it 
was found that at the onset of turbulence the resonance 
frequency would increase and from their figure we obtain 
v c i=s 62 mm/s. 2. A citation by the Osaka group [1(3] 
of their result with a vibrating wire. (This group has 
recently obtained more data [ll|.) 3. One result on vi- 
brating grids w 50 mm/s (l2| . We have no information 
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FIG. 2: (color online) Critical velocity for the onset of turbu- 
lence in superfluid 4 He as a function of the frequency. Blue 
squares: larger sphere; red circles: smaller sphere, tempera- 
ture is below 1 K; red triangle: vibrating wire at 7 mK [j|; 
blue inverted triangle: vibrating wire at mK temperatures 
[3~ot | ; x : vibrating grid at 8 mK [12|; diamond: tuning fork at 
1.5 K [l]|; + : tuning fork at 1.3 K [3]. The dashed line is 
a fit to all data (except for the +): v c = 2.0 (k2tt/) - 46 . 



on the accuracy of those data. It should be mentioned 
that because the velocity amplitude of these structures 
varies from zero at the fixed ends to a maximum in the 
middle, the measured velocity signal is some average from 
which the maximum velocity has to be calculated by some 
model. The resulting data are fairly consistent with ours. 
3. One result with a tuning fork of the Helsinki group 
[IH and one from Prague [3]. The calibration of the 
velocity scale at the tip of a tuning fork has been pre- 
sented in Ref . [l5| . We note that the value from Prague 
is lower than the result from Helsinki by a factor of ~3.5. 
The origin of this discrepancy is not clear, we can only 
speculate that signal losses in the cables might have oc- 
curred and remained undetected. Therefore, at present 
this data point remains uncertain to us and is not in- 
cluded in the fit. Fitting a power law to all the other 
data gives v c — 2.0(kw) - 46 (in SI units). It is obvi- 
ous that this fit gives convincing evidence of the scaling 
v c ~ yjnuj that was inferred before from purely dimen- 
sional considerations. If we fix the power to 0.5 we get 
as the final result v c ss 2.6 ^kuj. It should be mentioned 
that a very recent work by the Osaka group on wires vi- 
brating at various frequencies also shows the same scaling 
of v c , including the frequency range from 1 kHz to 10 kHz 



where we have no data 11|. 



QUALITATIVE MODEL 

The dimensional considerations leading to the fre- 
quency dependence of v c are supported by the follow- 
ing qualitative but very general argument. A "superfluid 
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Reynolds number" Re s = vI/k has been denned, where 
I is some characteristic length scale 16|,|l7|. F° r R- e s >1 
turbulence is expected to be stable when mutual friction 
is negligible. It is reasonable to assume that for oscilla- 
tion amplitudes a = v/lu that are small compared to the 
size of the sphere, a (< 30 fim) is the characteristic length 
scale. Hence, we may expect to observe turbulence for 
v > v c ~ \/kw ■ 

Another qualitative theoretical analysis is based on 
comparing the time scale of the oscillatory flow uj~ 1 with 
the time scale governing the dynamics of the vortex tan- 
gle. In a recent work on the onset of superfiuid turbulence 
in steady counterflow by Kopnin 17[ a rate equation for 
the length of vortex lines per unit volume L (m -2 ) has 
been derived: 



L = f3[v s L 3/2 ~nL 2 



(1) 



where — A(l — a') — Ba is determined by the coef- 
ficients of mutual friction a' and a and the constants 
A, B ~ 1 (more precise values for A and B would require 
numerical simulations with some particular geometry), 
and v s is the velocity of the superfiuid (with respect to 
the quasiparticle gas). In superfiuid 4 He ~ 1 except 
very near the Lambda transition. Equation (|T|) is a gen- 
eralization of the well known Vinen equation [l8| that 
determines the balance between the rate of growth of L 
due to vortex multiplication and the losses due to mu- 
tual friction with the quasiparticles. Stable solutions of 
Eq. (JTJ) are L—0 and the saturated value L sat = (v s /k) 2 . 
An analytical solution for > and constant v s can be 
found in . A simple approximation is derived by Kop- 
nin for the limit of an initial growth from a small value 
L < L sat at t=0: 



L/L sat = [{L sat /L Q ) 1/2 - («/t)]- 



(2) 



where the time constant r = 2k//3v 2 sets the time scale 
for the initial growth. After a time t ~ t*J L sa tJ La a very 
rapid growth ("vortex avalanche") sets in. Finally, the 
stable saturated value L sat is approached exponentially: 
L sat - L ~ exp(-i/r). 

In order to extend this scenario to the case of oscilla- 
tory flow we first mention the assumptions made in de- 
riving Eq. ([T]) and discuss whether these assumptions will 
hold in our case as well. Firstly, it is assumed that the 
normal fluid is at rest with respect to the reference frame, 
i.e., in our case with respect to the sphere. This is qual- 
itatively fulfilled because the normal fluid is roughly at 
rest with respect to the sphere within the viscous pen- 
etration depth that is of the same order as the radius 
at about 1 K while the superfiuid moves with respect to 
the sphere. Secondly, the superflow is assumed to be ini- 
tially homogeneous. The superfiuid velocity field around 
the sphere varies over a distance of the order of the ra- 
dius R while the vortex spacing is given by L~ x / 2 <C R, 
when values of L sat at typical v c are used [20(. Therefore, 



this condition is satisfied as well. Finally, some remarks 
concerning the initial value of Lq are in order. Suppose 
the drive is increased from zero without any remanent 
vorticity in the cell (Lq = 0), then L will remain zero, as 
indicated very impressively by the recent Osaka measure- 
Typically, in the experiments it is difficult to 
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ments 

prepare a helium sample that is completely free of re- 
manent vortices, because even natural radioactivity and 
cosmic rays are known to produce them Therefore, 
ultimately some strongly trapped vortex loops will be un- 
pinned and L will grow. In addition, the oscillating body 
may shed new vortices at sharp protuberances where the 
local flow velocity may be greatly enhanced. All that 
may occur at velocities larger than v c and that is why 
during an initial up sweep of the driving force, v c can be 
largely exceeded as shown in Fig. [TJ On the other hand, 
when reducing the drive in the turbulent regime, L is 
large and Lq is irrelevant. Below 0.5 K, where the switch- 
ing between turbulent phases and potential flow occurs 
when approaching v c , we assume that enough vorticity 
will have survived during the laminar phases, either on 
the surface of the body or somewhere else in the cell, so 
that the potential flow will break down again statistically 
and a new turbulent phase is born. This will continue un- 
til too little power from the drive is available to sustain 
a turbulent phase at all and potential flow will be stable. 

The central result of our work is the frequency depen- 
dence of v c . This we analyze now on the basis of Eq. ([1]) 
and Eq. ((2|) as follows. The time scale r determines 
how fast L can reach the avalanche where the very rapid 
growth sets in. If L is to reach the avalanche, this must 
happen on a time scale shorter than one period. Quali- 
tatively, we take the time in which the velocity changes 
from zero to the maximum, i.e., about one quarter of a 
period. Hence, we have the condition 



UJT < 1/4, 



and therefore 



> v c w ^8kuj/0 = 2.83^kuj/0. 



(3) 



(4) 



To compare Eq. ^ with the experimental data we eval- 
uate for 4 He by referring to tabulated values of a and 
a' Setting the constants A, B = 1, we find = 0.95 
(at 1.3 K), 0.89 (at 1.6 K), and 0.79 (at 1.9 K). Below 1 
K mutual friction is very small, hence we set 0=1. This 
implies a slow increase of v c by about 10% in qualita- 
tive agreement with experimental results as displayed in 
Fig. El 

Note that the source driving the oscillatory flow, in 
particular the geometry, the size, and the surface prop- 
erties of the oscillating body, do not appear in Eq. Q, 
i.e., we have a true universal intrinsic critical velocity 
below which turbulence cannot exist in oscillatory flow. 

We have also investigated numerically the evolution 
of the vortex density in oscillatory flow. A simple first 
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FIG. 3: (color online) Temperature dependence of the crit- 
ical velocity of a 100 /im sphere oscillating at 236 Hz. The 
variation due to (3(T) is indicated by the dashed line. 



choice might be to let v s in Eq. (HJ vary harmoni- 
cally. This, however, periodically changes the sign of the 
growth term which obviously is incorrect. Alternatively, 
one could try v s (t) ~ |sin(cjt)|, where the modulus is 
used to keep the correct sign of the growth term. But 
this is a unidirectional motion and not an oscillation. It 
is only useful to confirm that the initial value L > 
is irrelevant once L has begun to grow and that L(t) 
can follow the time dependent velocity field v s (t) only if 
condition is fulfilled, otherwise the oscillations of L 
are much reduced and L sat cannot be reached any more. 
Obviously, this approach is too simple and a rigorous ex- 
tension of Eq. {U to the oscillating case is needed but, 
to our knowledge, is presently not available. 



SIMULATIONS 

We have also simulated the oscillations of a sphere in 
the superfluid at zero temperature using the vortex fila- 
ment model. Assuming a smooth sphere, the calculations 
with a R=100 /im sphere indicate a critical velocity of 
order 100 mm/s at 200 Hz (see also Ref. Q), which 
is about 2. . .3 times higher than the experimental one. 
Those time demanding simulations are still far from be- 
ing complete. Preliminary calculations indicate that the 
critical velocity goes up with frequency. Also the calcu- 
lations for temperature 1.5 K, where the normal flow is 
taken as Stokes flow, predict that the temperature de- 
pendence of the critical velocity is weak, like estimated 
from the (3 term. It is clear that a single vortex ring 
in homogeneous oscillating flow is stable. Therefore the 
boundaries and inhomogeneous flow around the object 
are the main reasons why the initial (remanent) vortices 
go unstable. But this does not mean that there cannot 
be a "vortex avalanche" later when the line density has 
grown large enough containing enough vortices of differ- 
ent size, which is the situation on which Kopnin's model 
is based. It is also true that the simulations indicate 
that the initial increase in vortex density occurs during 
several periods, but same occurs also with Eq. ([T]) if ini- 



tially Lq <SS L sat . Initially the increase in vortex length 
is due to vortex stretching and subsequent collision to 
the sphere, like noted in Ref. [2l|]. Unfortunately, we 
have not managed to reach a steady state due to limited 
computer resources available. The range of scales that 
can be simulated is also very limited. It is still an open 
question what is the turbulent state really like around 
the vibrating objects, whether it is a state where the vor- 
tex density is of order L sat — (v s /n) 2 or much smaller. 
For that, one should take into account the sphere motion 
self-consistently, like done by Kivotides et al. 23j. In 
principle near v c the line density should be low, but so 
far the simulations only indicate either a fast grow (con- 
sidered as v > v c ) or drop to some small line density 
(with no clear velocity dependence) if we start from con- 
figuration containing few hundred vortices. Hysteresis is 
also seen in simulations. If initially the vortex number 
is small, and especially with small vortices, the required 
velocity, above which the turbulence develops, is larger 
than if started with the configuration with many vortices. 



CONCLUSIONS 

In summary, our results can be understood, at least 
in a qualitative manner, both by employing the super- 
fluid Reynolds number and by applying Kopnin's theory 
to oscillating superflow. A more rigorous theoretical de- 
scription of the onset of turbulence in oscillatory flow is 
desirable. Hopefully, also further progress in the numeri- 
cal simulations will improve our understanding of the ex- 
perimental results. Finally, we note that Eq. ^ can be 
applied to turbulence in superfluid 3 He as well, provided 
the temperatures is low enough so that (3 > can be 
achieved 24} and the oscillation frequency is low enough 
so that v c is smaller than the Landau critical velocity for 
pair breaking. 
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